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We study the spin-wave dispersion of localized spins in a disordered double-exchange 
model using the perturbation theory with respect to the strength of the disorder potential. 
We calculate the dispersion upto the next-leading order, and extensively examine the case of 
one-dimension. We show that in that case, disorder yields anomalous gapped-like behavior 



c/2 ' at the Fermi wavenumber of the conduction electrons. 
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Q '_ Colossal magnetoresistance manganites have been providing hot topics owing to their 

rich magnetic and electronic properties.^"^ Recent studies have clarified that the double- 

^ . exchange (DE) mechanism alone is not enough to understand a variety of phases of these 

^.f-^ ■ materials, which is actually induced by the interplay between spin, charge and orbital degrees 

^_. _ of freedom. 

Cp I Recently, the softening of the spin-wave dispersion, ^^"^^ observed in Ri-^^D^^MnOs (R=La, 

'NT , 

O ■ Pr, Nd and D=Ba, Sr, Ca, etc.), has suggested that a random potential also plays a crucial 

j^ ■ role in magnetic properties of these materials. ^^"^^ Although some other mechanism^^"'^^ have 

been proposed to describe the softening phenomena, recent numerical calculations by Motome 

'w ' and Furukawa^^"^^ have shown that a simple DE model with a random potential can explain 

O ■ the anomalous behavior of the spin-wave spectrum including broadening, gap formation as 



well as softening. 

Motivated by their work, we study a disordered DE model by the use of the perturbation 



H ' with respect to the strength of a random potential. Based on a general formula we derived, 

we show that in the case of 1-dimensional (ID) systems randomness yields a singularity of 
the spin-wave dispersion at the Fermi wavenumber of conduction electrons. It turns out that 
such anomalous behavior is mainly due to the nesting of the Fermi surface in ID models. We 
then discuss the higher dimensional systems, based on the results of ID systems. 
We start with the following Hamiltonian: 

where s, = ^cJ^o-o-o-'Qcr' and Si denotes, respectively, the electron spin and the localized spin in 
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the spin-S" representation, and e^ is the quenched random potential with statistical properties 
ej = and ejey = gdij. This model includes two large parameters, S and Jh- Assuming that 
localized spins are almost aligned to the z-direction, we first take into account the leading 
terms with respect to the series of 1/S. To this end, we utilize the Holstein-Primakoff mapping 
for localized spins: 

S+ = V2Sai, Sr = V2Sal, S^ = S - alai, (2) 

which leads us to the following truncated Hamiltonian up to 0{1/S^): 

H = -t Yl 4aCja - -^ Yl ^4aCia 
(i,j),a i,(T 



Y - J// Y'^44^Cil + tticl^Ci^) + ^Y1 4('i'^4a<^i<r + Y ' 



J 



(3) 

t,cr 

Based on this Hamiltonian, we calculate the spin-wave dispersion perturbatively with respect 
to the strength of disorder g, collecting the leading terms with respect to Jh- 

Previously, Furukawa^'^ calculated the magnon dispersion in the clean DE model, according 
to the diagrams (a) in Fig.l. In the presence of disorder, these diagrams should be evaluated by 
the use of the electron propagator with the self-energy induced by the scattering of electrons 
by impurities. Since down-spin electrons are located in the momentum space far above the 
Fermi energy in the large Jh limit, it is readily seen that disorder yields the self-energy to 
up-spin electrons only, which is evaluated as 

St(ie„) = Y -■ — I — ~ -sgn(e„) (4) 

in a weak disorder limit, where e„ = (2n -|- 1)tt/P is the fermionic Matsubara frequency, 
Efc = —2t ^ cos kfj, — £f is the bare electron dispersion, and the relaxation time r is given by 
1/r = 7rgD{eF)- 



(a) 





(b) ( = -^ — ' < - < + < '<■<"<-< + ... 

Fig. 1. (a) The leading self-energy of the spin- wave. Thin-lines denotes the down-spin propagator, 
whereas the thick-line denotes the up-spin propagator with the self-energy induced by disorder, 
(b) The propagator of up-spin electrons used in the diagrams (a). The thin- line and dotted- line 
denotes the bare up-spin propagator and impurity potential, respectively. 
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Using the propagator with the self-energy above, the diagrams in Fig.l (a) are calculated as 



Ilq (itOn) 



Jr^ \^ v^ 

2 ^^ lEr, 



-1 



-1 



Jh 



EE 



-1 



2 ^-^ "-rf iSr. 






25 



where 



1 / 2 

/^ (e) = — I 1 arctan re 

2 V vr 



(5) 



(6) 



and the arrow in eq.(5) implies the large Jh limit. This equation indicates that disorder merely 
smears the Fermi surface at the leading order and hence, does not give rise to anomalous 
behavior for the magnon spectrum. 

In order to calculate the self-energy of the spin-wave at the next leading order, let us 
consider the diagrams in Fig. 2. Each diagram is evaluated as follows: 



iP 

q+p q 





(C) 



Fig. 2. Diagrams which contribute to the next leading order of the spin- wave self-energy. The dotted- 
line denotes the impurity potential. The thick-wavy-line stands for the spin-wave propagator in 
Fig.l. 



r(l«) 



nq+p,q(^w„;ep) 



4^^pEE 



-1 



-1 



-j^ i£n - £k + S|(ien) i^n - £k+p + S^(ier 



i{£n + UJn) — £k+p+q - JrS 

2iu>n + 2efc — £k+p+q ~ ^k-q 



-^^pXp+TT^^P^fri^k 



£fc ~ £fc^ 



(7) 
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Sn k 



i{en + Wn) — efc+p+q - JhS 

]^^pY.fr{eul (8) 



2^ . 



and 



Jh 



2 ^ ^ ie„ - Efc + S|(ie„) ze„ - Sk+p + St(^^n) 



2 ^pXp, (9) 

where ep stands for the impurity potential in eq.(3) in the momentum representation, Xp is 
defined by 

^^ = Y^ fr{ek)-Jr{ek+p) ^ (10) 

and ujn is the bosonic Matsubara frequency a;„ = 2n7r//3. Collecting these diagrams, we obtain 

n?|p,,(^-n;ep) = ^Y. friek) '^"" + "'- " "^-'^ + "'-^^ " "^^^^'^ . (11) 

As already claimed, the contractions of the impurity lines in the same electron loop in 
Fig. 2, which corresponds to Fig.l, do not give any anomalies in the magnon spectrum. There- 
fore, one can expect singular behavior in another type of diagrams, in which impurity lines 
are contracted among different electron loops in such a way as induces the correlations among 
up-spin electrons through the disorder potential. One of such simplest diagrams is given in 
Fig. 3. To calculate this diagram explicitly, we use the quadratic dispersion e^ ~ tk'^ — ep. 

P 

...--X-.-..N 

q q+p q 

Fig. 3. The next-leading self-energy of the spin-wave. 

In ID systems, this approximation works basically without loss of generality, whereas in the 
case of higher dimensional systems, it is well-known that the above self-energy, the function 
Xp in particular, is quite sensible to the form of the dispersion relation of bare electrons. 
Therefore, we restrict our calculations to ID systems to reveal the mechanism of gap opening 
in a disordered DE model. As to higher dimensional systems, we will discuss what we expect 
at the end of the calculation based on the results of the ID model. 
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The use of the quadratic dispersion leads to the following expression of the self-energy in 
Fig. 3 for the ID system, 



n 



(1), 



lUJr. 



p 



{pq + q^) xl 



lUJr. 



n 



(0) 



(12) 



Here, we have neglected the Matsubara frequency in the numerator, since it merely gives rise 
to small effects. The total self-energy is now given by Ilg (iu;„)+ng (iu;„), and the spin- wave 



dispersion relation can be determined by solving lo =XVq (a;) + IT, 
solve this nonlinear equation iteratively. 

At the first step, one finds the spectrum given by 

nw = 2t'(i 



«, 



UJ] 



In what follows, we 



Wr, 



COS (7), 



(13) 



where t' = t/{2S) ^^ fri^k) cos k. Next, substituting this into 11, 



(1)/ 



[UJ 



we have 



LOn 



n(o) ^ ^ \^ {pQ + <ffxl 
" + 52 Z^ (0) _ j,(o) • 



n(o) _ ^ v [(^'?)' + '^'] 4 

<? SH' ^ ' 

p 



p^ + 2pq 



(14) 



r(0) 



~ t'q\ 



where we have used the quadratic dispersion approximation for the spin-wave, Ilg 

The second term above should exhibit singular behavior in the spin-wave spectrum. To 
observe this, we first calculate the function Xp with the smeared distribution function fri^) 
in eq.(6). We obtain 

{If +pkF cos ^ + kl 



1 1 1 

Xv = In 

^^ 27rt 2p 



(i)^ 



pkp cos ^ + kp 



(15) 



where kp and 9 is defined, respectively, kp = {kp + e^)^''^ with e = l/(iT) and tan^ = e/kp. 
It may be instructive to take the zero disorder limit, 5 ^ (e ^ 0) in this function: It then 
becomes the well-known form, 



1 1, 
Xp^ -77-:- In 



p + 2kF 



(16) 



p — 2kF 

This function plays a crucial role in the spin-density-wave (SDW) in ID systems due to the 
divergence at p = 2kF in the absence of disorder. However, in the system under consideration, 
it shows just a peak atp = 2kp due to the self-energy of electrons, i.e, the scattering of up-spin 
electrons by the disorder potential. In Fig. 4, we show examples of this function in two kinds 
of the hole densities (of the lower band) x defined by /ci? = 7r(l — x). One can see that the 
hight of the peak depends weakly on x, since g enters Xp through 9. 

Having derived the function Xp, we next calculate the dispersion of the spin-wave. The 
p-integration in eq.(14) is rewritten as 



E 

p 



{p^q^ + g^) xl 
p^ + 2pq 



1 



(27rt)2 47r 



^'Q' 



(2) ^ Q' 



.(0) 



(17) 
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0.5 1 1.5 2 2.5 3 



p/2kp 

Fig. 4. The function x in cq.(15) computed for two kinds of holc-dcnsitics: x ~ 0.5 and x = 0.3. The 
parameters used are g/t^ = 0.1 and 5* ~ 3/2 in both cases. 



where Q = q/kp and 

(n) _ J_ 

9q - 23 



An) 



for n = 0,2 is defined by 
1 1 \ 1 



dx 



+ 



r3—n 



In^ 



1 + 2 cos I + x^ 



2 cos I + x2 



^x + Q x-Qj x-^-"- \ 1 
Using the contour integral in the complex plane, this integral can be converted into 



(18) 



(n) 

9q 



vrwcos 



Q' 



^5, 



+7r 



X Ln 



n,0 



dxlm. 



cos ■ 




(x + ism^)^ " 



2i sin 2 + 2^ 



cos ■ 



2i sin ■ 



2 -""J^^^ 2 ■"/ 

This expression is convenient for numerical computation as well as to take the g 
as we shall see momentarily. The dispersion of the spin-wave is finally given by 



UJ„ 



t'q^ 



9 



1 kp 



(2) 

9q 



4 



, • (0) 
+ — 5q 



(19) 
limit, 

(20) 



An) 



SH' (27r)2 47r 

Computing Qq^ in eq.(19) numerically, we show in Fig. 5 examples of the spin-wave dispersion 
for the same systems as those in Fig. 4. One can indeed find an anomalous dispersion at the 
Fermi wavenumber of electrons. Although gap opening is not observed, a sharp jump appears 
just at the Fermi wavenumber. Therefore, we can conclude that the diagram in Fig. 3, which 
includes the correlations among up-spin fermions through disorder, plays an important role 
in the dynamics of localized spins, even though the disorder is nonmagnetic. 

So far we have derived the spin-wave dispersion in the presence of disorder by the use 
of /r(e) in the calculation of Xv ^^^ hence, of gn ■ To clarify the origin of such anomalous 
behavior, let us take g ^ limit in eq.(19), as has been done for Xp in 6q.(16), yielding 



(n) 

9q 



■K 



2vH 



In 



1 + 



for 




(21) 
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Fig. 5. The dispersion in cq.(20). The dotted-line denotes the quadratic dispersion in the 5 = hmit 
(i.e, ujq = t'q^). The down-arrows denote the positions of the Fermi wavenumbers. 

If one uses this formula instead of eq.(19), one reaches the fohowing spin- wave dispersion: 



UJn 



{t' + 2ga)q^ =F ^kpgaq In 



q + kp 



q-kp 



for 




(22) 



where the constant a is given by a = l/^GAirkpS'^t'). The dispersion in eq.(22) now diverges at 
q = kp. Therefore, the present result reveals that an enhancement in Xp at the nesting vectors 
is responsible for the anomalous gapped-like behavior of the spin-wave dispersion of localized 
spins. However, it should be stressed that even though the divergence in Xp is suppressed by 
disorder, which is due to the smearing of the Fermi surface, singular behavior of spin-wave 
dispersion seems robust. 

So far we have demonstrated that nonmagnetic impurities of conduction electrons actually 
induce the anomalous behavior of localized spins at the Fermi wavenumbers of electrons, 
as suggested by Motome and Furukawa, by using the leading order perturbation theory. In 
models with small Jh, the diagrams in Fig.l leads to interesting physics such as SDW or spin 
glass (RKKY exchange interaction). Contrary to these, the same diagram merely yields the 
canonical dispersion of the spin-wave in the present systems with large Jh because of weak 
correlations between electrons with different spins. Even in such systems, however, disorder 
induces correlations among up-spin electrons, and strongly affects the magnetic properties, 
even if disorder is weak and nonmagnetic. 

Finally, we would like to discuss the models in higher dimensions. The chief ingredient in a 
step- function-like behavior in the dispersion is the function Xp, which has a sharp peak at p = 
2kF in ID. This reflects the nesting of the Fermi surface. In higher dimensions, however, this 
function is continuous if we use the free electron (quadratic) dispersion. Therefore, in such an 
approximation, one cannot expect any anomalous dispersions in the spin-wave excitation. On 
the other hand, the tight-binding electron dispersion gives rise to an enhancement around the 
nesting wavevectors of electrons near half-filling. Therefore, one can expect similar anomalous 
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behavior of spin-waves there even in higher dimensions. The question is how robust such 
anomalous behavior is away from half-filhng. Unfortunately, cosine dispersion of electrons 
does not allow analytic calculations in higher dimensions, and therefore, the details, including 
numerical results, will be published elsewhere. 

The authors would like to thank J. Igarashi, Y. Motome and M. Yokoyama for valuable 
discussions. This work was supported in part by Grant-in-Aid for Scientific Research from 
JSPS. 
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